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Homogeneous neutron matter at subnuclear densities becomes unstable 
towards the formation of inhomogeneities. Depending on the average value of 
the neutron density one can observe the appearance of either bubbles, rods, 
tubes or plates embeded in a neutron gas. We estimate the quantum correc- 
tions to the ground state energy (which could be termed either shell correction 
or Casimir energy) of such phases of neutron matter. The calculations are 
performed by evaluating the contribution of the shortest periodic orbits in 
the Gutzwiller trace formula for the density of states. The magnitude of the 
quantum corrections to the ground state energy of neutron matter are of the 
same order as the energy differences between various phases. 
PACS numbers: 21.10.Dr, 21.65.+f, 97.60. Jd, 05.45.Mt 

It is somewhat surprising to find out that the problem we are going to present here, which 
could have naturally emerged as a typical study case in quantum chaos, has somehow eluded 
the attention of this community. Various billiard problems, in particular the Sinai billiard 
in 2D, are paradigmatic in quantum chaos studies. The structures believed to appear in the 
crust of neutron stars, namely bubbles, tubes, rods and plates, known also as the "pasta 
phase" 0], are nothing else but one of several natural realizations of Sinai billiards. The 
energetics of these and similar phases at higher densities involving quarks has been studied 
for many years |l| (this reference list is not meant to be exhaustive). 



All the studies mentioned above have estimated the ground state energy of neutron 
matter only in the liquid drop or Thomas-Fermi approximation. In this approximation 
applied to finite fermi systems - such as an atom, a nucleus, a metallic cluster or a quantum 
dot - one fails to put in evidence the fact that in the case of closed shells these systems show 
an enhanced stability. Apparently, an agreement has been reached in literature concerning 
the existence of the following chain of phase changes as the density is increasing: nuclei 
— > rods — > plates — > tubes — > bubbles — > uniform matter. There are a couple of studies, 
where the role of quantum corrections has been taken into account, one being the Hartree- 
Fock calculations of Refs. M. In Ref. || the shell effects due to the bound nucleons in 
nuclei immersed in a neutron gas only (mainly protons) have been computed. There it was 
concluded that this type of quantum corrections to the ground state energy will not lead 
to any qualitative changes in the sequence of the nuclear shape transitions in the neutron 
star crust. However, no other type of quantum corrections was suggested. We have shown 
that the quantum corrections to the ground state energy of neutron matter arising from 
the unbounded motion of fermions outside the bubbles in particular, are more important 
quantitatively and thus can lead to qualitative changes in the physics of neutron matter at 
subnuclear densities |3J]. In particular, since the magnitude of the quantum corrections to 
the ground state energy is of the same order of magnitude as the energy differences between 
various "pasta phases", there is no reason to expect that the sequence "nuclei — > rods — > 
plates — > tubes — > bubbles — > uniform matter" is occurring and instead strong disorder could 
dominate over various regular lattices. 

Let us start our analysis with reviewing a somewhat simpler case, that of systems with 
only one bubble |5||| . About half a century ago it was predicted that if a large nucleus was 
ever to be observed, most likely it would have either the shape of a torus or of a sphere with 
an empty spherical cavity inside (this is what is typically referred to as a bubble nucleus) [|7| . 
Such geometries lead to a minimum of the sum of Coulomb and surface energies. Recently 
it was suggested that many fermion systems with a similar geometry can be created in the 



case of highly charged atomic clusters |§. In order to focus our attention on the quantum 
corrections, we shall discard here the role of long range Coulomb interaction, which is easy 
to account for, and we shall sidestep the issue of bubble stability as well and refer the 
interested reader to earlier literature. In the case of atomic clusters the Coulomb energy 
plays a secondary role anyway || . For the sake of simplicity, we shall consider only one type 
of fermions with no electric charge. 

It is known that in the case of a saturating iV-fermion system the total energy has the 
general structure 

E(N) = e v N + e s N 2 '* + e.N 1 / 3 + E SC (N). (1) 

The first three terms represent the smooth, liquid drop part of the total energy. E SC (N) is 
the pure quantum contribution, known as shell correction energy, the amplitude of which 
grows in magnitude approximately as oc iV 1 / 6 , see Ref. Q. The liquid drop part of the 
energy in the case of billiard systems is referred to as the Weyl energy p!0| , pT| . 

One of the simplest questions one can raise is: How should one determine the position 
of a bubble or void? Considerations based on naive symmetry arguments have been used 
in the past, without any critical consideration, to simply position the bubble in the very 
center of the system. Recently we have shown that moving the bubble off-center can often 
lead to a greater stability of the system, due to shell correction energy effects alone . It 
is worth noting at this point that the shell correction energy bears a remarkable similarity 
to the Casimir energy in quantum field theory and in critical phenomena [ 12- 14| . The 



Casimir energy can be thought of as a specific measure of the magnitude and nature of 
the fluctuations induced in the energy spectrum by the presence of various "obstacles" and 
is computed with a formula very similar to the one used to estimate the shell correction 
energy ||. There is a large number of studies of various fermion systems, which shows that 
the character of the shell corrections is strongly correlated with the existence of regular or 



chaotic motion in the classical limit [Tl,iq-tt8|. Integrable systems have conspicuous large 



gaps in their spectra, while chaotic systems have a spectrum characterized by rather small 



fluctuations and essentially no large gaps. If a spherical bubble appears in a spherical system 
and if the bubble is positioned at the center, then for certain magic fermion numbers the 
shell correction energy E SC (N), and hence the total energy E(N), has a very deep minimum. 
However, if the number of particles is not magic, in order to become more stable the system 
will in general tend to deform. Real deformations lead to an increased surface area and 
thus to an increased liquid drop energy. In the case of systems with bubbles/voids there is 
a somewhat unexpected extra soft deformation mode - bubble positioning. Merely shifting 
a bubble off-center deforms neither the bubble nor the external surface and therefore, the 
liquid drop part of the total energy of the system remains unchanged. Therefore, one would 
expect that classically this deformation mode has exactly zero energy. Quantum mechanics 
at first glance seems to lend some support to this idea. In recent years it was shown that 
in a 2-dimensional annular billiard, which is the 2-dimensional analog of spherical bubble 



nuclei, the motion becomes more chaotic as the bubble is moved further from the center [19 



One might anticipate then that the importance of the shell corrections diminishes when the 
bubble is off-center and thus becomes easier to move a bubble, the further from the center it 
is. One might also naively assume that the most favorable arrangement is the one with the 
bubble at the very center, when the system has the highest symmetry. It came as somewhat 
of a surprise to find out that this is not the case most of the time and that the profile of the 
energy surface has a very unexpected structure |5|U . 

We shall refer the reader to our previous work for more details concerning the physics 
of systems like bubble nuclei and revert now to the case of bubbles, rods, tubes and plates 
floating in an infinite neutron gas. All calculations we present here for infinite systems were 
performed at constant chemical potential. In order to better appreciate the nature of the 
problem we present here, let us consider the following situation. Let us imagine that two 
spherical identical bubbles have been formed in an otherwise homogeneous neutron matter. 
One can then ask the following apparently innocuous question: "What determines the most 
energetically favorable arrangement of the two bubbles?" According to a liquid drop model 



approach the energy of the system should be insensitive to the relative positioning of the 
two bubbles. A similar question was raised in condensed matter studies, concerning the 
interaction between two impurities in an electron gas. In the case of two "weak" and point- 
like impurities the dependence of the energy of the system as a function of the relative 
distance between the two impurities a is given by 

E(a) = ljdrij rfr 2 V 1 (r 1 ) X (r 1 - r 2 - a)V 2 (r 2 ), (2) 

where xi r i ~ r 2 — a ) is the Lindhard response function of a homogeneous Fermi gas and 
Vi(ri) and V2(r 2 ) are the potentials describing the interaction between impurities and the 
surrounding electron gas. At large distances, when k F a 3> 1, this expression leads to an 



interaction first derived by Ruderman and Kittel |20|j21 



. cos(2k F a) . . 

£(a) cx 1 / \ (3) 

where k F is the Fermi wave vector and m is the fermion mass. This asymptotic behavior 
is valid only for point-like impurities, when k F R <C 1, where R stands for the radius of 
the two impurities. This condition is typically violated for nuclei and bubbles immersed in 
a neutron gas, for which k F R 3> 1. We have shown in Ref. Q that in the case of large 
"impurities" (k F R S> 1) and at constant chemical potential the interaction energy at large 
separations a becomes 

2m V a J 7r 2 

The most striking aspect of this interaction is the fact that it is inversely proportional to the 
square of the separation, instead of 1/a 3 as in Eq. (|]). The bubble-bubble interaction has 
a surprinsingly long range. Apart from Coulomb and gravitational interactions, there are 
no other longer ranged interactions known between two objects. One can advance a simple 
qualitative explanation of the difference in the power law exponent of the impurity-impurity 
interaction (0) and the bubble-bubble interaction in terms of waves. A large "obstacle" 
reflects back more of the incident wave than a point object, which acts like a pure s-wave 



scatterer. We shall see below that perfectly planar surfaces lead to an interaction which 
decays with distance even slower, inversely proportional to the separation between surfaces. 
Naturally, when the incident and reflected waves interfere constructively, one does expect 
a more energetically favorable geometric configuration. Obviously, the mere fact that the 
interaction (|j) oscillates shows that for each value of a = ir(n + 3/4) /kp, where n is an 
integer, the potential has a minimum and one thus expects that a bubble-bubble molecule 
with such a radius can be formed. Moreover, for each such size various vibrational and 
rotational states, corresponding to small fluctuations of the bubble-bubble distance and 
angular velocity are most likely to exist. Tunneling between various molecular sizes should 
also appear. However, the treatment of such molecular systems and tunneling should be 
deferred until the bubble inertia is also estimated. 

The arguments of the cosine in Eq. (0) and of the sine in Eq. @ are nothing else but 
the classical action in units of % of the bouncing periodic orbit between the two impurities 
and this suggests the most likely method of attack for this problem, the semiclassical ap- 
proximation based on the Gutzwiller trace formula for the density of states. As a matter of 
fact this is the way Eq. (f|) has been derived ||. The formation of various inhomogeneities 
in an otherwise uniform Fermi gas can be characterized by several natural dimensionless 
parameters: kpa 3> 1, where as above a is a characteristic separation distance between two 
such inhomogeneities; kpR ^> 1, where R is a characteristic size of such an inhomogeneity; 
and kps = 0(1), where s is a typical "skin" thickness of such objects. The fact that the first 
two parameters, kpa and kpR, are both very large makes the adoption of the semiclassical 
approach natural. Since the third parameter, kps, is never too large or too small, one might 
be tempted to discard the semiclassical treatment in the case of saturating fermion systems 
altogether as unreasonable. However, there is a large body of evidence pointing towards the 
fact that even though this last parameter in real systems is of order unity, the seemingly 
unreasonable approximation kps <C 1, which we shall adopt in this work, is surprisingly ac- 



curate [11]. The corrections arising from considering kps = 0(1) lead mainly to an overall 



energy shift, which is largely independent of the separation among various objects embedded 
in a Fermi gas. On one hand, this type of shift can be accounted for in principle in a suitably 
implemented liquid drop model or Thomas-Fermi approximation. On the other hand, the 
semiclassical corrections to the ground state energy arising from the relative arrangement 
of various inhomogeneities have to be computed separately, as they have a different physi- 
cal nature. We are thus lead to the natural assumption that a simple hard-wall potential 
model for various types of inhomogeneities appearing in a neutron fermi gas is a reasonable 
starting point to estimate quantum corrections to the ground state energy, see Refs. flj6|.|TT 



and earlier references therein. One might expect that such simplifications will result in an 
overestimation of the magnitude of quantum corrections to the ground state energy, but 
that the qualitative effect would be reproduced. 

One can compute the quantum corrections to the ground state energy analytically for 
plate-like nuclei or voids with the neutron gas filling the space between slabs. The quantum 
correction to the energy for such a system is given by: 
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where the exact and the Weyl (smooth) energy [|K],|11|] per unit volume are given by 
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stands for the integer part of the argument in the square brackets, and a = L — 2R is the 
distance between slabs and R is half of the width of the slab. Here L 3 is the volume of an 
elementary (cubic) cell and the factor "2" in front stands for the two spin states. Using 
these formulas, one can show that the shell correction energy behaves as 
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where G(x) is an approximately periodic function of its argument. For x > 1, G(x + 1) ~ 
with properties G(x — n) w — 1 and — 1 < G(^) < 0.5. Note that at small separations 
a the shell correction energy is attractive in character. 

In the case of rods or tubes and bubbles we resorted to the Gutzwiller trace formula to 
estimate the corrections to the density of states. We were interested in what in the nuclear 
physics lingo would be called the "gross shell structure" and we discarded the fine details. As 
this approach implies a certain amount of spectral averaging one needs to account only for 
the existence of the shortest periodic orbits. As we are dealing with an infinite system there 
are no discrete levels and correspondingly no shells. However, the presence of "obstacles" 
induces fluctuations in both spatial matter distribution and density of states, which thus 
lead to the appearance of what we term shell correction or Casimir energy. For the case of 
spherical voids there are 26 periodic orbits between the nearest neighbors of three different 
lengths 2Li = 2(L - 2R), 2L 2 = 2(L v / 2 - R) and 2L 3 = 2(L v / 3 - R) in the case of a simple 
cubic lattice. Thus the shell correction/Casimir energy and density are equal to: 



1 [V 3 
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The contribution due to one periodic orbit to the fluctuating part of the level density reads: 
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and hence we get: 
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where Ai = 6, A2 = 12, A3 = 8 respectively. Above, the summation over n is over repetitions 
of a given orbit and such summation becomes superfluous when one performs a spectral 
averaging in order to extract the "gross shell structure" . This shell energy can be regarded 
as the interaction energy between bubbles immersed in neutron gas. The bubble-bubble 
interaction (Q) was obtained from the above formula for the case of two bubbles only in the 
asymptotic limit. One can derive similar formulas for the case of rods Q. The fact that the 
interaction energy between various homogeneities is described in terms of periodic orbits, 
points to another nontrivial aspect: the interaction between three or more such objects 
cannot be reduced to simple pairwise interactions. Orbits which bounce among three or 
more objects would lead to three-body, four-body and so forth genuine interactions. Such 
orbits were not yet included in our analysis, even though they can be accounted for in a 
straight forward manner. Since periodic orbits of this kind are typically longer, one expects 
their contribution to be masked when one performs a spectral averaging in order to extract 
the "gross shell structure". 

Lack of space prevents us from presenting more detailed results and discussing many 
other aspects of the energetics of inhomogeneous fermi systems. Hopefully, the few results 
we have highlighted here will convince the reader of the richness of these systems and of 
the spectacular role played by geometry and the chaotic versus integrable character of the 
single-particle dynamics. Besides static properties, one should expect that the dynamics of 
such systems will be extremely rich and challenging to describe. 



[1] G. Baym et al, Nucl. Phys. A175, 225 (1971); D.G. Ravenhall et al, Phys. Rev. Lett. 50, 
2066 (1983); M. Hashimoto et al, Prog. Theor. Phys. 71, 320 (1984); K. Oyamatsu et al, 



Prog. Theor. Phys. 72, 373 (1984); J.M. Lattimer et al, Nucl. Phys. A432, 646 (1985); R.D. 
Wilson and S.E. Koonin, Nucl. Phys. A435, 844 (1985); M. Lassaut et al, Astron. Astrophys. 
183, L3 (1987); K. Oyamatsu, Nucl. Phys. A561, 431 (1993); CP. Lorenz et al, Phys. Rev. 
Lett. 70, 379 (1993); C.J. Pethick and D.G. Ravenhall, Annu. Rev. Nucl. Part. Sci. 45, 429 
(1995); G. Watanabe et al, see Los Alamos e-print archive astro-ph /000 1 273 ; H. Heiselberg 
et al, Phys. Rev. Lett. 70, 1355 (1992). 

[2] J.W. Negele and D. Vautherin, Nucl. Phys. A207, 298 (1973); P. Bonche and D. Vautherin, 
Nucl. Phys. A372, 496 (1981); Astron. Astrophys. 112, 268 (1982). 

[3] K. Oyamatsu, M. Yamada, Nucl. Phys. A578, 181 (1994). 



[4] A. Bulgac and P. Magierski, Los Alamos e-print archive astro-ph/0002377 



[5] A. Bulgac et al, in Proc. Intern. Work, on Collective excitations in Fermi and Bose systems, 
eds. C.A. Bertulani and M.S. Hussein (World Scientific, Singapore 1999), pp 44-61 and Los 
Alamos e-preprint archive, nucl-th/9811028] . 



[6] Y. Yu et al, Phys. Rev. Lett. 84, 412 (2000). 

[7] H.A. Wilson, Phys. Rev. 69 538 (1946); J. A. Wheeler, unpublished notes; P.J. Siemens and 
H.A. Bethe, Phys. Rev. Lett. 18, 704 (1967); CY. Wong, Ann. Phys. 77, 279 (1973); W.J. 
Swiatecki, Physica Scripta 28, 349 (1983); W.D. Myers and W.J. Swiatecki, Nucl. Phys. A 
601, 141 (1996). 

[8] K. Pomorski and K. Dietrich, Eur. Journ. Phys. D 4, 353 (1998). 

[9] V.M. Strutinsky and A. G. Magner, Sov. J. Part. Nucl. Phys. 7, 138 (1976). 

[10] R.T. Waechter, Proc. Camb. Phil. Soc. 72, 439 (1972); H.P. Baltes and E.R. Hilf, Spectra 
of Finite Systems, Wissenschaftsverlag, Mannheim, Wien, Zurich: Bibliographisches Institut, 
(1976). 

[11] M. Brack and R.K. Bhaduri, Semiclassical Physics, Addison-Wesley, Reading, MA (1997). 



[12] H.B.G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948); V.M. Mostepanenko and N.N. 
Trunov, Sov. Phys. Usp. 31, 965 (1988) and references therein. 

[13] M.E. Fisher and P.G. de Gennes, C.R. Acad. Sci. Ser. B 287, 207 (1978); A. Hanke et al, 
Phys. Rev. Lett. 81, 1885 (1998) and references therein. 

[14] M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 (1999) and references therein. 

[15] O. Bohigas et al, Phys. Rev. Lett. 52, 1 (1984). 

[16] R. Balian and C. Bloch, Ann. Phys. 67, 229 (1972). 

[17] H. Nishioka et al, Phys. Rev. B 42, 9377 (1990); J. Pedersen et al, Nature 353, 733 (1991); 
M. Brack, Rev. Mod. Phys. 65, 677 (1993) and references therein. 

[18] V.M. Strutinsky, Sov. J. Nucl. Phys. 3, 449 (1966); Nucl. Phys. A 95, 420 (1967); ibid A 122, 
1 (1968); M. Brack et al, Rev. Mod. Phys. 44, 320 (1972). 

[19] O. Bohigas et al, Phys. Rep. 223, 43 (1993); O. Bohigas et al, Nucl. Phys. A 560, 197 (1993); 
S. Tomsovic and D. Ullmo, Phys. Rev. E 50, 145 (1994); S.D. Frischat and E. Doron, Phys. 
Rev. E 57, 1421 (1998). 

[20] M.A. Ruderman, C. Kittel, Phys. Rev. 96, 99 (1954). 

[21] A.L. Fetter and J.D. Walecka, Quantum Theory of Many Particle Systems, (McGraw-Hill, 
New York, 1971). 



